On some mathematical model of turbulent flow with intensive 

selfmixing 



M. BURNAT 

University of Warsaw 
Institute of Applied Mathematics and Mechanics 

O ' Banacha 2, 02-097 Warszaw, Poland 

^ ■ 

The paper presents a new possibihty of mathematical modeling of turbulent flow 

r^ I with the intensive selfmixing. The flow is described in this model by a family of 

mappings 

[^. (^t,^ -uj ^ S{t,f,u); u,S C R^,t>t° 

I I such that in general 

^ : cji n w^ = ^ S{t, t°, co^) n S{t, t", w) = 

The model allows a new approach to certain problems of turbulent flow. For ex- 
ample there are possibilities to introduce some models of turbulent chemical reactors, 

^ '. or description of the unpredictable explosions of turbulence in the calm laminate flow. 

\Q [ In this model we introduce mathematical tools that allow deflnition of various mixing 

"^ I states of the fluid. It is possible to formulate in the appropriate mathematical language, 

the idea that the given fluid can have various mixing states. It is also possible to de- 

Q ■ termine the influence of a given mixing state on the flow parameters such as velocity 

fleld and density. The first formulation of the model is given in [1] . 



% 



l.The general feature of the model 

For the balls 

oosix) = {y e R^ -.{y-xl <e},e> 0, 

consider the sets Ve(x,t) C -R^ of velocities of fluid particles moving at the time t in 
Us{x). The velocity sets may be measured. The laminar model and our model may be 
shortly characterized in the following way. 

Laminar model 

1) The flow is described by family of regular invertible mappings (diffeomorphisms): 

^^^o : u ^^ S{t, t", u); uS{t, t°, u) C i^^ t> t°. 



where S{t, t°, u) is the set containing at t > t° the whole fluid contained at t° in u, and 
no another one. It means that the fluid contained in u at the time t does not mix with 
any fluid in the neighbourhood, and keep for ever its identity. 

2) There exist the Euler velocity fleld v{x,t), 

(1) v{x, t) = lime^o Ve(x, t). 
Our model . 

I. The mappings 

^t,to ■.cu^S{t,t°,uj) 

describing the flow are not invertible. The set S{t,t°,u) is the minimal set containing 
at the time t > t° the whole fluid m{u,t°), but it may contain another fluid portions. 

II. Instead of the fleld of Euler velocities (1) we assume only the existence of fleld 
of velocity sets 

V(x,t) Ci^^ 

such that for numbers A > 0, £ > small enough the following approximation 

(2) S{t + A,t,Ueix))^x + AV{x,t) 

holds in the sense of measure. 

Basing on the assumptions I, II we can formulate the integral conservation laws of 
mass, impulse, energy and impulse momentum. Moreover, we obtain the equivalent 
closed integro-differential system. In this way we obtain a model of fluid flow in which 
different fluid portions are mixing one with another loosing their identity. 

To this end we have to introduce flrst the notion of the physical a-quantities. We 
begin with the deflnition of the a-density. This is a nonnegative function 



Q-.R'^R, g= g{x,t,a) ^ ■. ^ ^ d3 



> 0, aeY{x,t), 
0, aeR^\V{x,t) 



with the following properties. For the mass portion fllling up u at the time t, the 
following equalities hold 



\m{uj,t)\= I g{x,t)dx = K, / / g{x,t,a)dxda, 

Jlj Jw{x,t) 



Q{x,t) = /^ / Q{x,t,a)da 

Jw{x,t) 

where Q{x,t) is the usual density, and /t > is some scaling constant. 



We shall construct an approximation of g{x, t, a) explaining its physical sense. 
The sets 

S{t + A, t, UJe{x)) = $t+A,t(w), 

where A > 0, e > small enough, are the minimal sets that contain at the time t + A 
the whole fluid fllling up ujs{x) at the time t. This fluid has in S{t + A,t,Us{x)) at 
t + A density 

g{y) for y e S{t + A, t, Ueix). 

Using (2), we have 

\m{ue{x),t)\= / g{y)dy^ / g{y)dy = A^ g{x + aA)da. 

J S(t+A,t,uie{x) Jx+AV(x,t) Jv(x,t) 

Hence 

Q{x,t)^ ' \^\^' ^' ^— - / g{x + aA)da, 

and we obtain the following approximation of the a-density 

, ^ , ( g{x + Aa), aeV{x,t) , 3 /A\3 

^ ' ' ^ \ 0, ae R^\V{x,t) ' ATiKeJ 

In our model we use the general notion of the a - quantity. For example we deflne the 
a - impulse as ag{x, t, a), so that the impulse of the fluid portion m{u, t) is equal to 



imp{uj,t) = K, / ag{x,t,a)dxda 
Jw Jw{x,t) 

If some physical expression contains the k-fold integration over the variable a, then we 

multiply the integral by h!'. 

Given the a-quantities we may obtain the observed mean Euler velocity v{x,t) of 

the flow: 

imp{u,t) Jj^3ag{x,t,a)da 
v{x,tj = lim — 



d(^)^o \m{uj,t)\ f^^ g{x,t,a)d. 



a 



In almost all physical situations one may take some ball A G R^, for which we know 
that g{x, t,a) = for a G R^\A, and we may integrate over A only. 

2. The Mass Conservation Law 

We have to introduce the fluid portions m{u,a,t),a C A, which are parts of the 
fluid in u moving at the time t with the velocities a E a G A. We have 



■m{uj,a,t)\ = K, / g{x,t,a)dxda, \m{u, A,t)\ = \m{uj,t)\. 




Moreover, we introduce the mass mixer: 



M:R^^^R, M = M{x,t,a,P), 
which describes the following mixing processes: 



m(u,a,t) ^^ miuj,h,t) 



^'^ L la Ib\M(^^ ^' "' P)dxdad/3 



(3) 



m(u,a,t) ^^ m(oj,b,t) 

^ ' ' ^ a,bcA — ^ ' ' ^ 



K^ f^ /^ / V ^ M{x, t, a, /3)dxdad/3 



If in the mixing process one of the portions m{u, a, t),m{u, b, t) is underlined, then it 
means that we ask amount of mass (positive or negative equal to the integrals in the 
frames on the right) which is transported to the underlined portion in the unit time . 
It turns out that 

M(x, t, a, (3) = -M{x, t, (3, a) 

The following three mixing processes (4), (5), (6) should be taken into account in the 
mass conservation low for the portion m{u,a,t), where a G A: 



(4) 
(5) 
(6) 



m{u, a, t) <^=^ m{uj, A\a, t) 
m{u, a, t) <^=^ m{B?\(jj, a, t) 
m{uj, a, t) <^=^ m{R^\oj, A\a, t). 



The amount of mass transported in the unit time in this processes to the portion 
m(a;, a, t) is respectively equal to 



(4a) 



(5a) 



(6a) 



K^ f^ f^ f^.M{x, t, a, /3)dxdad(3 



— /t/Q^/^[< a,n{x)g{x,t,a) > —E < n(a;),grad^f)(x, t, a) >]dxda 



"'^^ hu, la lA\a < ^(^)' ^(^' ^' "' 1^) > dxdadjS 



where n{x) is the external normal unit vector to the boundary du. 

In (4) we used the properties (3) of the mass mixer M{x,t,a, (3). We consider the 
process (5) like in the laminar model without any mixer, so as in this process we have 
on the both sides of the boundary du, fluids with the same kinematical characteristics. 
In the process (5) we take into account the mass transport of the diffusion type with 
the constant E > 0. The process (6) is a mixing process of two fluid portions with 
different kinematical characteristics across the boundary du. Therefore we have to 
introduce the boundary mass mixer B, 



B-.R'"^ R\ 



B = B{x,t,a,/3). 



It describes the amount of mass transported in the unit time to m{uj, a, t) in the process 
(6). 

Finally taking together (4a), (5a), (6a) we obtain the following integral mass con- 
servation law for the portion m(cu, a, t). 



(7) 



LLY'tQix^t, 


a)+ < a 


grad^ 


g> 


—EAxg{x,t, a)+ 


'^JA\a 


divxB{x 


,t,a,P)- 


-M{x 


t,a 


/3) df3\dxda = 










u C R^,a C A 



We shall use the following simple 



Localization Theorem 



Let D : R^ -^ R, F : R^ -^ R. Then the condition 



dxda = 0, 



(8) fJ^[D{x,t,a) + Kf^^^F{x,t,a,/3)df3 

u C R^,a G A 
is equivalent to the following system of equations: 

(9) D{x,t,a) + Kj^F{x,t,a,l3)dl3 = 0, 

X e R^,a e A 



(10) 



F{x, t, a, /3) + F{x, t, /3, a) = 0, 

X e R*,a,P e A 



Applying the Localization Theorem to integral conservation law (7), we obtain the 
following integro-differential system equivalent to the integral law: 



;ii) 



dtQ{x, t, a)+ < a, gmd^g > -EA^g + n J^ 


dw^B{x, t, a, (5) - M{x, t, a, /3)j d(3 = 


X e R^,ae A 


div^B{x, t, a, (3) - M{x, t, a, /3) + div^B{x, t, /3, a) - M{x, t, /3,a) =0 


X e R^,a,/3 e A 



3. The Impulse Conservation Law 

We introduce the following impulse mixer J{x,t,a, P). It is the function 

J : R^^ ^ R^, J{x, t, a, /3) = aM{x, t, a, /3) + i{x, t, a, /3) 

where aM{x, t, a, (3) describes the transport of impulse caused by the mass transport, 
and i{x,t,a, P) describes other types of impulse changes. The physical sense of the 
mixer J{x,t,a, /3) is the following. The mixer J{x,t,a, (3) describes the amount of 
impulse transported in a unit time to the portion m{u, a, t) in the process 

^ ' ' ^ a,bcA V 5 5 7 

This amount is equal to 

(12) K J^ J^J^^^J{x,t,a, I3)dxdadl3. 



On the one hand, the change of the impulse in a unit time in the portion 772(0;, a, t) 
is equal to 

ndt / / ag{x,t,a)dxda. 



On the other hand, this change is equal to the force acting on m{u, a, t). This force 
is the sum of the following forces 



The forces acting on m{u, a, t) by fluid portions: 

m{uj, A\a, t) , m{R^\u , a, t),m{R^\u, A\a, t). 



+ 



The external forces acting 
on m{uj^ a,t). 



Let us determine this forces. The forces acting on m{u, a, t) are equal to the amount 
of impulse transported in a unit time to m{(jj, a, t). 



In the process 

m(a;, a, t) <^==^ m{uj, A\a, t) 

(according to the definition of the impulse mixer) the amount of impulse transported 
to m{uj^ a, t) in a unit time is equal to 

1^^ I I I J{x,t,a, j3)dxdadl3. 

Jul J a J A\a 

In the process 

m{(jj, a, t) <^==^ m{R^\u, a, t), 
the amount of impulse transported in a unit time to m{u, a, t), is equal to: 

dxda 



-K I I a < a, n{x) > g{x, t,a) ~ E < n{x), grad^g > 

' dui J a 

In the process 

m(co', a, t) <^=^ m(R^\a, A\a, t) 
the increase of impulse in the unit time in the portion m{uj^a, t) is equal to 

—K^ Iff JBix,t,a, l3) ■ n{x)dxdad/3 



doj J a J A\a 



where we introduced the boundary impulse mixer JB{x,t,a, f3). It is a matrix Jb '■ 
i?^° — 7> i?^ X i?^. Finally, let us take into account the processes (4), (5), (6), and the 
external forces. Then, after changing the surface integrals into volume ones, we obtain 
the following integral form of the impulse conservation law for the portion m{u, a, t): 




dtag{x, t,a) + a < a, grad^g > —EA^g + 
(13) +K, J^,[dwa:JB{x, t, «, /3) — J{x, t, «, (5)]d[5 — f{x, t, a) >dxda = 0, 

u C R^,a C A 

where 

F(t) = K, / f{x,t,a)dxda 

J ui J a 

is the external force acting on m(w, a, t). In the case of the earth gravitation field, we 
have 

f{x,t,a) = gg{x,t,a), (? = const. 



Applying to (13) the Localization Theorem, we obtain the following integro-difFerential 
system equivalent to the integral impulse conservation law: 



(14) 







dtag{x 


,t,a) + Q 


[< a,grad,j.Q 


> -EA^g = 
X e R^,ae A 






= 


kJ^[J{x, 


t,a,/3)- 


d[VxJB{x,t,0 


',/3)]dp + f{x,t,a) 




J{x 


t,a 


/3) + J{x 


t,/3,a)- 


-div^[JB(x,t, 


a,f3) + JB{x,t,f3,a)] = 
X e R^;a,/3 e 



A 



4. The Energy Conservation Law. 

Let us introduce the a-inner energy of the flow. This is a function 

e : R' ^ R,e = e{x, t, a), 

e{x, t.a) = for a G A\V{x, t). 

Denote by e{u, a, t) the inner energy of the fluid portion m{uj, a, t). Then the physical 
sense of the a-energy is explained by the formula: 



e(co',a, t) = K, / e{x,t,a)g{x,t,a)dxda. 




Denoting by e{u, t) the inner energy of the fluid portion m{u, t), we obtain the following 
equality: 

e{uj,t) = / e{x,t)g{x,t)dx = e{uj, A,t) = K, / / e{x,t,a)g{x,t,a)dxda. 

J LO J U) J A 



Hence, we get 

(15) 



e{x^ t)g{x, t) = K j^ 5(x, t, a)g{x^ t, a)da. 



The force F{u, a, t, ) acting on the portion m{u, a, t) is equal to the increase of the 
impulse of m{u,a,t) in the unit time. Hence, from the impulse conservation law we 
obtain 

F{uj,a,t) = K, / E{x,t,a)dxda, 

Jul J a 

where 

E{x,t,a) = K, [J{x,t,a)—diVxJBix,t,a,l3)]d/3+a[— < a, grad^g > +EAxg]+f{x,t,a) 

J A\a 



The increase of the energy of the portion m{u, a, t) in the unit time is equal to the 
power due to the force F{u,a,t). Let us evaluate it. Divide u and a into small parts 
Au C uj,Aa C a, and notice that F{u,a,t) acts on m{uj,a^t) is such a way that on 
the small parts m(Aw, Aa,t) C m{u,a,t) the following forces are acting: 



F{Auj, Aa,t) = K / / E{x,t,a)dxda. 

J Auj JAa 

The power due to F{u,a,t) is the sum of the powers P{Auj,Aa,t) due to the forces 

F{Au,Aa,t) where Au C u,Aa C a. For xe Auj,aE Aa, we obtain the following 
approximation for small At: 

P{Auj,Aa,t) ^ -r— ( I (^ I At) ( /t|Acj||Aa| )|i?(a;, t,Q;)| cos ^ {a,E{x,t,a)) = 



O „ / o o , 



= /t|Aa;||Aa| <a, E{x,t,a) > 
Finally, for the whole power P{u, a, t) due to the force F{u, a, t), we obtain 

P{uj,a,t) = y^ P{Auj, Aa,t) = k, < a, E{x,t,a) > dxda. 

Now, we can write down the integral energy conservation law for the portion m{uj, a, t): 



1 2 e{^,t,a) 
2 



< a, E{x, t, a) > >dxda = 0, 



(16) X; /a{^< [^(^' ^' ")^(a^> t^ a) + \a\ 

to C R^,aC A 
Putting a = A and taking into account (15), we obtain the following relation 



(17) dt{e{x,t)g{x,t)) = K f^ < a, f > —\a\'^l< a,gradxg > —EAxg + dt 



(){x,ta) 



da. 



In this way, one determine the mean inner energy e{x, t) for given a-density g{x, t, a). 
Applying to (16) the Localization Theorem, we obtain the following integro-differential 
system equivalent to the integral energy conservation law: 



a 



e{x, t, a)g{x, t, a) + \a 



2e 

2 



+ \a\ 



< a^grad^g > —EA^g = 

X e R^,a e A 



f^ < a, J{x, t, a, (3) — div^ 7^(3;, t, a, (5) > df3+ < a, f{x, t,a > 

< a, J{x, t, a, (3) — divj, Jb(x, t, a, /3) > + 
+ < /3, J(x, t, /3, a) — divj: 7^(2;, t, /3, a) >= 



5. The Impulse Momentum Conservation Law. 



Applying the principle 



The derivative dtoi the impulse momentum 



The momentum of acting forces 



and our Impulse Conservation Law we can formulate the Impulse Momentum Conservation 
Law for the fluid portion ?t2(w, a, t). 

If we denote FAH = {{FAH)i, {FAH)2, {FAH)3), for F,H e i?^ then we obtain 
the following integro-differential system equivalent to the integral Impulse Momentum 
Conservation Law 



(19) 



dt{x A ag)j = div^, E{x A a)jgradxQ — a{x A a)j —{x A /(x, t, cy))j+ 



+fi Jj^ [x A J{x, t, a, (5) j +dxi {x A h^)j + d^^ {x A h'^)j + d^.^ {x A h^)j d(3 

j = l,2,3,x e R^,ae A 
(x A J{x,t,a,(3)] +dj:^(xAh^{x,t,a,(3)] +dx^{x A h'^)j + d3:.^{x A h^)j + 

+ (x A J{x,t,(3,a)j + dxi(x A h^{x,t, (3,a)] +(92^2(3; A h'^)j + dx.^{x A h^)j = 

j = l,2,3,a:G R^,a,f3 e A 



where the matrix Jb is denoted as 



Jb{x, t, a, /3) = (h)), h\x, t, a, /3) = {h\, h^, h\) 



6. The constitutive relation and the full closed system of the model. 

In order to formulate the constitutive relation let us express approximately the 
boundary mixer B{x,t,a, (3) as a function of the mixer M{x,t,a, /3). To this end, 
introduce the following notations 



d'oj = [j \^y:\y 



x\ < e 



D^'u = cfuj\uj C R\uj, (fu = cfu n w C w, 
where e is a positive small number. 
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Consider two following complex mixing processes (1°) and (2°). These processes 
taken together, give us an approximation of the process (6): 

m{uj, a, t) <^=^ m{Ft'\(jj, A\a, t) 



:i°) 



(a) 


m{D''u;,A\a,t) < 


J=^ m{D'^u, a,t) 


iP) 


■m{D'^u,a,t) <^=^ 


> m{(jj, a,t) 



The complex process 1° gives us: 

m{D^uj, A\a, t) <^==^ m{uj, a, t). 

This process represents a part of the process (6). The second complex process is of the 
form: 



(2°) 



(a) 


m{d'^u},a,t) <^=^ 


> m{d'^u, A\a,t) 


(b) 


m{d^u,A\a,t) < 


^^^m{R^\uj,A\a,t) 



This complex process gives us 

m{d^ijj^ a, t) <^==^ m(R\uj, A\a, t). 

Hence the complex processes (1"), (2°), represent together an approximation of the 
process 

771(0;, a, t) <^=^ m{R\uj, A\a, t) 

let us analyse more precisely the processes (1°) and (2°). 

The process 1°. 

The amount of the mass transported in the unit time in the process (1°), (a) to 
m{D'^u, a, t) is equal to 

K^ / / M{x,t,a, I3)dxdadl3 

JD^iuj) J a J A\a 

Hence the density q{x, t) and the a-density q^x^ t, a) of this mass are equal to 

g{x,t) = K^ / / M{x,t,a, /3)dad/3, 

J a J A\a 



11 



and 



Q{x,t,a) = K, / M(x,t,a, /3)d/3. 

lA\a 



The process 1°, (b) is of the type (5). Replacing in the description of the process 
(5) £>(a:,t, a) by k, J^,M{x,t,a, /3)df3 and neglecting for simplicity the mass transport 
of diffusion type, we obtain that the amount of mass transported in the process (1°), 
(b) in a unit time to m{uj, a, t) is approximately equal to 

—K^ / <n{x),a> / M(x,t,a, l3)d/3 dxda 

J dui J A\a 

2 



= —K / / / < n{x),aM{x,t,a, l3) > dxdad/S. 

J du) J a J A\a 

The process (2°). 

The density of mass transported in the unit time in the process (2°) (a) to m{d'^u, A\a, t) 
is equal to 

g{x,t) = K^ / / M{x,t,a, I3)dadl3, x G d'^u, 



J A\a J a 

and the a-density of this mass is equal to 

q{x, t,a) = k / M(x, t, a, (5)d(5, x G d^to, a G A\a. 

J a 

The process (2°)(b) is of the type (5), but in the opposite direction from oj to R^\u. 
Hence neglecting for simplicity the mass transport of diffusion type the amount of mass 
transported in the unit time in the process (2°)(b) to the portion m{R^\u,A\a,t) is 
approximately equal to 

K^ / <n{x),a> M{x,t,a, I3)dl3 dxda = 

J du) J A\a J a 

= K^ j I I < n{x), l3M(x,t,a, l3) > dxdad/3. 

Juj J a J A\a 

In the final description of the mass transport to m{uj,a,t) in the process (6) we 
have to take the last expression with sign minus because it describes the mass learving 
tu. Considering together the complex processes (1°) and (2°), we see that the amount 
of mass transported in the unit time in the process (6) to m{u,a,t) is approximately 
equal to 

—K^ / / / < n{x),aM{x,t,a, l3) + l3M{x,t, l3,a) > dxdad/S. 

J dui J a J A\a 
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Hence we obtain the following approximated expression for the boundary mixer: 

(20) B{x, t, a, 13) ^ aM{x, t, a, l3) + l3M{x, t, l3, a). 

We have the integro differential system of 16 equations: (11) 2 equations (the mass 
conservation law), (14) 6 equations (the impulse conservation law), (18) 2 equations 
(the energy conservation law), (19) 6 equations (the impulsemomentum conservation 
law). The system contains 18 unknown functions, a-quantities, and mixers: 

g, e, M, Bi, B2, B3, Ji, J2, Js, Jb^; i,j = 1,2,3. 

Basing on the relation (20), we introduce the following constitutive relation : 

(21) B{x,t,a,l3)= aM{x,t,a,l3) + l3M{x,t,l3,a) b{x,t,a,l3), 

where 

b:R^^-^R, b = b{x,t,a,l3), 

is a new unknown function. In this way we reduce the number of the unknown functions 
to 16, and we obtain the general closed system of the model. 

7. Simplified Models. 

The full integro-differential system of the model is very complicated and we are 
forced to seek for some simplified models. There exist many possibilities of such sim- 
plification. We shall give here the simplest of them. The main idea is to suggest some 
special form of the mass mixer M{x, t, a, /3), for example 

(22) M{x, t, a, 13) = $(«, /3)/i(a;, t, a, /3), 
where $(q;, (3) > 0, $(q;, /3) = $(/3, a), 

^{a, 13) = cos-, i9=^{a,l3), 



fi(x,t, a, (3) = 
D = const > 0, and 



g{x,t, a)r{Dd) for d = \/3\g{x,t, /3) — \a\g{x,t, a) > 0, 
g{x, t, /3)r{Dd) for d = \(3\g{x, t, (3) — \a\g{x, t,a) <0 



, -f-. for d>0 
^(^) = < ^^' for d < 0. 



i-d 

A short justification of this form of mixer M{x,t,a, f3) is as follows. Consider 
for a, P & A two small fluid portions approximately moving with the velocities a 
and f3. The mixer M{x,t,a, P) describes the amount of mass (positive or negative) 
transported in the unit time from /3 to a. Hence for the mixer (22), if li = |/3|^(x, t, (3) — 
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\a\g{x,y,a) > then the portion moving with velocity a loses (in the unit time), in 
favour of the portion moving with the velocity /3, the amount of mass proportional to 

cos —g{x, t, a)r{Dd). 

\i d = |(/3)|f)(a;, t, /3) — \a\Q{x^t^a) < 0, then the mass of a portion moving with the 
velocity /3 loses (in the unit time), in favour of the portion moving with the velocity 
a, the amount of mass proportional to 

cos —q{x, t, (3)r{Dd). 

Our simplest model is based on the following assumptions. 

1° We neglect the mixing process (6), m{u,a,t) <^=^ m{R^\ijj,A\a,t), putting 
B{x,t,a,(3) =0. 

2° We introduce the form (22) of the mass mixer M{x, t, a, /3), defining in this way 
the mixing state of the fluid. This is closely related to the impulse conservation law, 
so as our mixer (22) decides how the impulse situation governs the mixing process. 

3° In the simplifled model, since B[x,t,a, P) = 0, the mass conservation law takes 
the following form 



(23) 



4° We use the energy conservation law in the form (17). The solution (){x,t,a) of 
the closed system (23) and the thermodynamical state equation p = f{e, g) allow to 
construct the Euler parameters v{x,t), g{x,t),e{x,t),p{x,t). 

Notice, that one may consider many different forms of mass mixers M{x,t,a, P). 
For example, we may introduce some thermodynamical or electrodynamical parameters 
to the deflnition of M{x, t, a, (3) . 

Observe that for the flow described by the a-quantities the construction of the sets 

S{tX.^) = ^t,to{uj) 
must be based on the following relations: 



(9jf)(x, t,Q;)+ < a,gradx0 > 


-E/\^Q = kJ^ M{x, t, a, I3)d(3 
X e R^,a e A 



S{t + /\,T,UJs) 



|j[y + AV(|/,r) 



r > t°, ujs = {y : \y — x\ < 5}; hence 



the smaller A > 0, the better approximation. The value of A > should be chosen 
according to the properties of the sets V(x,t),t > that appear in the initial and 
boundary conditions. If V(a;,t) changes in x,t rapidly, then we have to take A > 
small enough otherwise A > may be larger. 
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The question of determination of parameters e > and k = J;(— )^ is as follows. 
First we must realize that without strict definition of k, the system (23) and the model 
are not defined in full. Taking k, = const larger or smaller we decide that the mixing 
step in the flow considered will be larger or smaller, respectively. 

For instance, this may be seen in our condition 

S{t + A, t, Ue) ^x + AV(x, t). 

The right side does not depend on e. Hence, for big k and small e the set u^ is very 
small, and the very small mass portion rniue, t) must be spread on the set a; + AV(a;, t) 
which does not depend on e. 

The first numerical experiments for our simplified model may be found in [2]. 
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